Prediction of Ferromagnetic Correlations in Coupled Double-Level Quantum Dots 
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Numerical results are presented for transport properties of two coupled double-level quantum 
dots. The results strongly suggest that under appropriate circumstances the dots develop a novel 
ferromagnetic correlation at quarter-filling (one electron per dot). In the strong coupling regime 
(Coulomb repulsion larger than electron hopping) and with the inter-dot tunneling larger than the 
tunneling to the leads, an S=l Kondo resonance develops in the density of states, leading to a peak 
in the conductance. A qualitative "phase diagram", incorporating the new FM phase, is presented. 
It is also shown that the conditions necessary for the ferromagnetic regime are less restrictive than 
naively believed, leading to its possible experimental observation in real quantum dots. 
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Strongly correlated electronic systems, such as high-T c 
cuprates, heavy fermions, and manganites, display a va- 
riety of nontrivial collective states, which are difficult to 
analyze due to the many-body character of the interac- 
tions, and the difficulties in experimentally controlling 
the parameters determining these interactions. These 
problems are severe in materials that spontaneously grow 
in particular structures and patterns, with several ef- 
fects (lattice, spin, charge, orbital) in direct competition. 
Therefore, the observation of a celebrated many-body ef- 
fect, the Kondo effect, in a single quantum dot (QD)i has 
captured the attention of the strongly correlated com- 
munity It is conceivable that the most interesting states 
that are spontaneously stabilized in some materials - and 
are very difficult to control - could instead be artificially 
created in a man-made structure. In this framework, a 
natural first-step is to analyze coupled QDs. In fact, the 
two-impurity Kondo problem - extensively studied since 
the 80's^ - can now be realized in a real system. More- 
over, recent investigations have reported antiferromag- 
netic (AF) correlations between two single-level coupled 
QDs, in competition with Kondo correlations^^. As a 
consequence, it is now clear that two of the most remark- 
able magnetic states known to exist in spontaneously 
grown materials - the Kondo and AF states - have already 
found realizations in the context of QDs. However, the 
other dominant magnetic state of some materials - the 
ferromagnetic (FM) state - has comparatively received 
much less attention^. For the dream of artificially repli- 
cating collective states using QDs to be fulfilled, a realiza- 
tion of FM states must be achieved. The lack of attention 
to FM states in QDs should not be surprising in view of 
the physics of FM materials, such as manganites. Here, 
the FM state is reached by removing electrons (doping) 
from an AF state. Under the constraint of having one 
particle per level (1/2-filling), and only one level active 
per QD as in most previous investigations, the double- 
exchangei generated FM state cannot be realized. To 
reach a FM state, more levels need to be active, result- 
ing in less than one electron per level. 



In this Letter, clear evidence is presented for the 
development of ferromagnetic correlations between two 
double-level QDs*: at 1/4-filling (one electron per dot), 
two coupled double-level QDs form a triplet state. Cou- 
pling this state to ideal metallic leads produces a Kondo 
resonance and a peak in the conductance. The results 
do not appear to be restricted to only a pair of QDs, 
but they seem valid for larger QD arrays. Basically here 
it is reported a realization of the double- exchange mech- 
anism using QDs. Although the above mentioned ef- 
fect is stronger if the appropriate intra-dot inter-level 
many-body interactions are added to the Hamiltonian, 
it is important to stress that these interactions are not 
necessary: considering just an intra-level Coulomb repul- 
sion (Hubbard U) is enough to obtain qualitatively the 
same results, opening the possibility for the FM regime 
to be experimentally observable. Figure 1 schematically 




FIG. 1: Schematic representation of the two double-level QDs 
system and the hopping amplitudes of Eqs. (1) and (3). 

depicts the system analyzed here and introduces the la- 
belling for the different tunneling parameters t' and t". 
To model this system, the impurity Anderson Hamilto- 
nian that describes the two QDs with two levels (denoted 
a and (3) is given by 
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where the first term represents the usual Coulomb re- 
pulsion between electrons in the same level (considered 
equal for both levels). The second term represents the 
Coulomb repulsion between electrons in different levels 
(the U' notation is borrowed from standard many-orbital 
studies in atomic physics). The third term represents the 
Hund coupling ( J > 0) that favors the alignment of spins 
and the fourth term is the energy of the states regulated 
by the gate voltage V g . To decrease the number of free 
parameters, all the calculations presented here assume 
the following relations: U' = 2U/3 and J = U —U'. As 
discussed later, the main result in this Letter does not de- 
pend on the specific values taken by U' and J. Note that 
a and /3 are separated by AV, and by modifying this 
parameter an interpolation between one- and two-level 
physics can be obtained. The last term represents the 
inter-dot coupling, with matrix element t" . For simplic- 
ity, we assume that there is no hopping between levels a 
and (3. The dots are connected to the leads (represented 
by semi-infinite chains) by a hopping term with ampli- 
tude t' , while t = 1 is the hopping amplitude in the leads 
(and energy scale). More specifically, 
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where c\ i(J (c\ irJ ) creates electrons at site i with spin a 
in the left (right) contact. Site "0" is the first site at 
the left of the left dot and at the right of the right dot, 
for each half-chain. The total Hamiltonian is Ht = + 
■fficads + ffint- Note that for V g =-U/2-U' + J/2-AV/2, 
the Hamiltonian is particlc-holc symmetric. Using the 
Keldysh formalism 9 , the conductance through this sys- 
tem can be written as 10 a = ^-\t 2 Gi r (E F )\ 2 [p^f)] 2 - In 
practice, a cluster containing the interacting dots and 
a few sites of the leads is solved exactly, the Green 
functions are calculated, and the leads are incorporated 
through a Dyson Equation embedding procedure (details 
of the embedding have been provided elsewher e^ 11 ). In 
Figs. 2-4, the cluster used involved the two QDs plus one 
lead site at left and right. 

In Fig. 2a, results for conductance (solid line) at strong 
inter-dot tunneling (t"/t' 3> 1) are shown. The main 
feature displayed is the peak at V g ~ —1.0 (at and near 
1/4-filling). The occupancy for this value of V g is ap- 
proximately 1 electron per dot (dashed line x4) and the 
total spin St of the four levelsi£(dotted line) is « 1.0. 
The smooth charging of the levels as the gate poten- 
tial decreases (in the peak region) indicates a possible 
Kondo regime. This is confirmed in Fig. 2b, where the 
density of states (DOS) close to the Fermi level is dis- 
played as the gate potential varies from 0.0 to —1.5 (top 
to bottom). Through this variation of V g the two dots 
are charged with one additional electron (the total mean 
charge varies from w 1.6 to 2.6). One can clearly see a 
Kondo resonance pinned to the Fermi level. For lower 
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FIG. 2: Results for the strong inter-dot tunneling regime 
it" /t' ^> 1) showing the presence of a ferromagnetic state at 
1/4-filling and concomitant Kondo effect: (a) Conductance 
(solid line), charge occupancy per level and per spin (dashed 
line), and total spin St in the four levels of the dots (dotted 
line) vs. gate potential V g for two coupled double-level QDs. 
The values for the parameters are indicated. Note that the 
maximum in the conductance corresponds to a spin 1 Kondo 
peak that occurs when there is exactly one electron per dot 
and the total spin St ~ 1.0. The other two peaks occur when 
there are 1 and 3 electrons in the two-dots system, corre- 
sponding to St = 1/2 Kondo resonances, (b) DOS indicating 
the presence of a Kondo resonance associated with the St = 1 
Kondo peak in (a). The numbers on the right side indicate 
the total number of electrons inside the two quantum dots. V g 
takes the values 0.0, —0.5, —0.83, —1.12, —1.5 from the top 
to the bottom curve. The Fermi level is located at ui/t = 0.0. 



values of the gate potential (in the region at and near 
1/2-filling, with 2 electrons per dot) the conductance is 
drastically reduced and the total spin St inside the dots 
reaches its minimum value, indicating the formation of 
a global singlet state. Calculations of the total spin in 
each dot indicate that this singlet state is formed by the 
AF coupling of two spins S = 1. A description of how 
this picture changes as the t" decreases is shown in Fig. 
3a, where results are shown for five different values of t". 
The conductance at the particle- hole symmetric point, 
V g = —5.0, varies from zero, for t" = 1.0, to 1.0 (in units 
of e 2 /h), for t" = 0.08. Fig. 3b shows how the Kondo cor- 
relation (between the total spin in the dots and a conduc- 
tion electron in the first site of one of the leads) evolves 
from a negligible value for t"/t' > 1 to a large value 
(w -0.65) for t"/t < 1. The inset of Fig. 3a displays 
the change of the total spin (from St ~ 0.0 to « 3/4) 
as t" decreases. The two main peaks in the conductance 
discussed up to now were the St = 1 Kondo peak at 1/4- 
filling (relevant in the strong inter-dot tunneling regime) 
and the peak at 1/2-filling (relevant in the weak inter-dot 
tunneling regime). It is interesting to discuss how these 
peaks evolve as AV increases. Fig. 4a shows results for 
t" jt' > 1 ft" = 1.0, U = 5.0 and t' = 0.2). The solid 
line displays the conductance and the doted line displays 
the total spin St- Level separation AV increases from 
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FIG. 3: (a) Conductance G vs. gate potential V g for different 
values of inter-dot tunneling (t = 1.0,0.6,0.4,0.2, and 0.08, 
see convention in (b)). V g varies across the region where the 
fourth electron is charged into the double-dot system. Note 
that the conductance at V g = —5.0 gradually increases from 
zero (at strong tunneling, t" = 1.0) to the maximum value 
(at weak tunneling, t" = 0.08). This variation indicates a 
transition from two S — 1 spins (in each dot) forming a 
global singlet to two uncorrelated S = 1 spins, each form- 
ing a Kondo resonance with the lead conduction electrons. 
The inset shows the variation with t" of the total spin inside 
the dots, (b) Variation with t" of the spin-spin correlation 
between the total spin in the double-dot and a conduction 
electron located in the first site of the leads. 



bottom to top (values for each graph are displayed in the 
left side). From AV = 0.0 up to AV « 0.6 the width of 
the conductance peak slowly decreases, as also does the 
maximum value of St- Above AV" ~ 0.7 (not shown) the 
narrowing of the peak accelerates (as does the decrease 
of St), until the peak has all but vanished for AV" = 1.0. 
For higher values of AV (top graph, AV = 10.0), the 
conductance shows the typical Coulomb blockade profile 
previously discussed 5 - for coupled single-level dots when 
t"/t' » 1. 

Figure 4b shows the corresponding results for t"/t' < 1 
(t" = 0.08, U = 5.0 and t' = 0.2). Note that the cen- 
tral peak does not change appreciably from AV = 0.0 to 
AV = 2.0. In fact, changes start only above AV = 3.0, 
when the central peak splits in two (AV > 3.2, not 
shown). For AV > 3.4 the two peaks start moving far- 
ther apart from each other and become very narrow. Fi- 
nally, for AV = 4.0 the central peaks have disappeared, 
and the remaining structures are already similar to the 
single-level result t"/t' < 1. The top graph (AV = 30.0) 
is basically the result reported for single-level QDs at 
weak interdot tunneling 5 (if one discards the slight shoul- 
ders in the internal peaks). 

Based on the results displayed on Fig. 4a, a qualitative 
phase diagram for the strong inter-dot tunneling regime 
can be sketched. In Fig. 5, the electron occupancy is 
in the horizontal axis (controlled by V^) and AV is in 
the vertical axis. The left side (indicating 1/2-filling) is 
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FIG. 4: Variation of the conductance with level separation 
AV. The level separation is indicated for each curve, (a) t" = 
1.0: The gate potential decreases down to the particle-hole 
symmetric value, to highlight how the St = 1 Kondo peak 
varies with A V. Note that all graphs have the same horizontal 
axis scale, except for the upper one (scale indicated on top), 
(b) t" = 0.08: The gate potential varies down to the lowest 
value (total charging of the dots), to highlight the variation 
of the central peak (at the particle- hole symmetric point). 
Here also the upper graph has a different horizontal scale. A 
discussion of how the data interpolate between double- and 
single-level dots is given in the text. 



dominated by antiferromagnetism for all values of AV < 
U. The singlet formed by the four levels is made of two 
spins S ~ 1.0. For AV > U one recovers the single- 
level picture. The right side of the phase diagram, which 
describes the evolution of the central peak in Fig. 4a, is 
more interesting. For AV < t" one has the novel FM 
region, where an St = 1 Kondo effect is present. For 
AV > t" an AF region is found, with no Kondo effect. 

A finite size scaling analysis was done (results not 
shown) to verify how our numerical results converge with 
cluster sizeii. The authors found out very little change 
in the results with increasing cluster size, giving us con- 
fidence that all the qualitative results here discussed are 
not caused by finite-size effects. It is also important to 
stress that the calculations presented in Fig. 2a were re- 
produced for U' = J = 0.0 (with the values of all other 
parameters kept the same as before (U = 5.0, t" = 1.0 
and t' = 0.2)) and the results obtained barely changed. 
This indicates that the FM correlation and the St = 1 
Kondo should be experimentally observable, since the 
only requirement is to have two double-level QDs with 
strong inter-dot tunneling^. 

Figure 6 qualitatively summarizes the main result pre- 
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FIG. 5: Qualitative phase diagram for the strong inter-dot 
tunneling regime (t" jt' ^> 1). Horizontal axis indicates the 
level occupancy (controlled by V g ) and vertical axis indicates 
the level separation (controlled by AV). In the 1/2-filling 
region, displayed in the left side, one goes from an AF coupling 
between two S = 1 spins (for AV < V g ) to a situation where 
the two lower levels are completely occupied (for AV > V g ). 
In both cases there is no Kondo effect. On the other hand, 
in the 1/4-filling region (right side), for AV < t" one has the 
novel St = 1 FM Kondo region, which gives way to an AF 
region (with no Kondo effect) once AV > t". The regime 
with 3 electrons in the two dots is very narrow as a function 



1 / 4- filling, FM correlations will develop and conductance 
through a Kondo channel is allowed, (b) On the other 
hand, single-level coupled QDs will develop AF correla- 
tions in the strong inter-dot tunneling regime and con- 
ductance is suppressed. The results discussed in this pa- 
per complete the analogy between QD states and mag- 
netic phenomena in bulk materials. Previous investiga- 
tions had shown that Kondo and AF states were possible 
in QDs. Now, at least theoretically, a regime with ferro- 
magnetism has also been found, if more than one level per 
dot is active. Certainly, it would be important to confirm 
experimentally this prediction. Our calculations empha- 
sizing multilevel dots present analogies with multi-orbital 
materials such as manganites, nickelates, cobaltites, and 
ruthenates. These compounds have a plethora of phases, 
all of which could find realizations in QDs systems as 
well. 
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Schematic representation of the main result in this 



sented in this Letter: (a) For double- level coupled quan- 
tum dots in the strong inter-dot tunneling regime at 
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Note that the total spin in the four levels in the two 
quantum dots (denoted as St and obtained through 

{St )ciustcr=5'T(5'T + 1), where St = X^.a^.*' wnere i 
labels the dots and A labels the levels) is not a good quan- 
tum number (since the dots are not isolated). However, St 
gives a good indication on the nature of the Kondo effect 
(S = 1/2 orS = 1). 

For a finite size scaling analysis of similar systems, see C. 
A. Biisser et al, to appear in Phys. Rev. B. 
One can understand why an effective FM J e f / is generated 
(even in the absence of the Hund J and the U' terms) by us- 
ing the following argument: At 1/4-filling (two electrons in 
the four levels), the hopping from one lead to the adjacent 
dot (and vice-versa) is maximized when the spins of the 
electrons in the levels a and /3 in the same QD are parallel 
to each other. A strong inter-dot tunneling then generates 
an effective FM coupling between the QDs. Therefore, the 
only restriction for the FM coupling seems to be that the 
leads should have just one conduction channel. 



